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ABSTRACT 

The present study, which is the first of a series of 
investigations dealing with specific issues of low 
pressure turbine (LPT) boundary layer aerodynamics, 
is aimed at providing detailed unsteady boundary flow 
information to understand the underlying physics of the 
inception, onset, and extent of the separation zone. A 
detailed experimental study on the behavior of the 
separation zone on the suction surface of a highly 
loaded LPT-blade under periodic unsteady wake flow 
is presented. 

Experimental investigations were performed at 
Texas A&M Turbomachinery Performance and Flow 
Research Laboratory using a large-scale unsteady 
turbine cascade research facility with an integrated 
wake generator and test section unit. To account for a 
high flow deflection of LPT-cascades at design and off- 
design operating points, the entire wake generator and 
test section unit including the traversing system is 
designed to allow a precise angle adjustment of the 
cascade relative to the incoming flow. This is done by 
a hydraulic platform, which simultaneously lifts and 
rotates the wake generator and test section unit. The 
unit is then attached to the tunnel exit nozzle with an 
angular accuracy of better than 0.05°, which is 
measured electronically. 


Utilizing a Reynolds number of 110,000 based on 
the blade suction surface length and the exit velocity, 
one steady and two different unsteady inlet flow 
conditions with the corresponding passing frequencies, 
wake velocities and turbulence intensities are 
investigated using hot-wire anemometry. In addition to 
the unsteady boundary layer measurements, blade 
surface pressure measurements were performed at Re= 
50,000, 75,000, 100,000, and 125,000 at one steady and 
two periodic unsteady inlet flow conditions. 

Detailed unsteady boundary layer measurement 
identifies the onset and extent of the separation zone as 
well as its behavior under unsteady wake flow. The 
results presented in ensemble-averaged and contour plot 
forms contribute to understanding the physics of the 
separation phenomenon under periodic unsteady wake 
flow. Several physical mechanisms are discussed. 

NOMENCLATURE 

c blade chord 

c ax axial chord p ~ p 

C p pressure coefficient, C = 

d R rod diameter ” ( Pi ~ P)inl 

M number of samples 

N number of wake cycles 

p, p, static, total pressure 

p, static pressure of first tap on suction surface 

Re Reynolds number Re = S 0 E 7 /v 
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S B blade spacing 

S R rod spacing 

s streamwise distance from blade leading edge 

s 0 blade suction surface arc length 
t time 

Tu turbulence intensity 

U belt translational velocity 

V ax axial velocity 

V exit exit velocity 

V velocity 

v fluctuation velocity 

y normal distance from the surface 

y blades stagger angle 

v kinematic viscosity 

•u one wake-passing period 

o cascade solidity, o = c/S B 

cp flow coefficient, cp = V ax /U g 

Q reduced frequency Q = — = — — 

s R v m <P S R 

INTRODUCTION 

In recent years gas turbine engine aerodynamicists 
have focused their attention on improving the efficiency 
and performance of the low pressure turbine (LPT) 
component. Research at industry, research centers, and 
academia has shown that reduction in blade count can 
be achieved without substantially sacrificing the 
efficiency of the LPT-blading. This reduction 
contributes to an increase in thrust/weight ratio, thus 
reducing the fuel consumption. Contrary to the high 
pressure turbine (HPT) that operates in a relatively high 
Reynolds number environment, the LPT experiences a 
variation in Reynolds number ranging from 50,000 to 
250,000 dependent on operation conditions. Since the 
major portion of the boundary layer, particularly along 
the suction surface, is laminar, the low Reynolds 
number in conjunction with the local adverse pressure 
gradient makes it susceptible to flow separation, thus 
increasing the complexity of the LPT- boundary layer 
aerodynamics. The periodic unsteady nature of the 
incoming flow associated with wakes substantially 
influences the boundary layer development including 
the onset and extent of the laminar separation and its 
turbulent re-attachment. Of particular relevance in 
context of LPT-aerodynamics is the interaction of the 
wake flow with the suction surface separation zone. 
While the phenomenon of the unsteady boundary layer 
development and transition in the absence of separation 
zones has been the subject of intensive research that has 
led to better understanding the transition phenomenon, 
comprehending the multiple effect of mutually 
interacting parameters on the LPT-boundary layer 
separation and their physics still requires more research. 


The significance of the unsteady flow effect on the 
efficiency and performance of compressor and turbine 
stages was recognized in the early seventies by several 
researchers. Fundamental studies by Pfeil and Herbst 
[1], Pfeil et al. [2], and Orth [3] studied and quantified 
the effect of unsteady wake flow on the boundary layer 
transition along flat plates. Schobeiri and his co- 
workers [4, 5, 6, 7] experimentally investigated the 
effects of periodic unsteady wake flow and pressure 
gradient on the boundary layer transition and heat 
transfer along the concave surface of a constant 
curvature plate. The measurements were systematically 
performed under different pressure gradients and 
unsteady wake frequencies using a squirrel cage type 
wake generator positioned upstream of the curved plate. 
Liu and Rodi [8] carried out boundary layer and heat 
transfer measurements on a turbine cascade, which was 
installed downstream of a squirrel cage type wake 
generator mentioned previously. 

Analyzing the velocity and the turbulence structure 
of the impinging wakes and their interaction with the 
boundary layer, Chakka and Schobeiri [7] developed an 
intermittency based unsteady boundary layer transition 
model. The analysis revealed a universal pattern for the 
relative intermittency function for all frequencies and 
pressure gradient investigated. However, the above 
investigations were not sufficient to draw any 
conclusion with regard to universal character of the 
relative intermittency function. Further detailed 
investigations of the unsteady boundary layer along a 
high Reynolds number turbine blade of the Space 
Shuttle Main Engine (SSME) by Schobeiri et al. [9] 
and its subsequent analysis [10], [11], and [12] verified 
the universal character of the relative intermittency 
function. For this purpose, Schobeiri et al. [9] utilized 
a conceptually different type wake generator, which is 
also used for the investigation presented in this paper. 
Fottner and his co-workers [13], [14] and Schulte and 
Hodson [15] used the same wake generating concept for 
investigating the influence of unsteady wake flow on 
the LPT-boundary layer. Kaszeta, Simon, and Ashpis 
[16] experimentally investigated the effect of unsteady 
wakes on laminar-turbulent transition within a channel 
with the side walls simulating the suction and pressure 
surfaces of a LPT-blade. They utilized a retractable 
cascade of cylindrical rods for generating the unsteady 
wakes. Lou and Hourmouziadis [17] experimentally 
investigated the effect of oscillating inlet flow 
conditions on laminar boundary layer separation along 
a flat plate under a strong negative pressure gradient 
which was imposed by the opposite wall. 
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Using the surface mounted hot film measurement 
technique, Fottner and his co-workers [13], [14] 
Schroder [18], and Haueisen, Schroder, and Hennecke 
[19] documented strong interaction between wakes and 
the suction surface separation zone on LPT- blades, 
both in wind tunnel cascade tests and in turbine rig. 
Furthermore, they investigated the boundary layer 
transition under the influence of periodic wakes along 
the LPT-surface and found that the interaction of 
wakes with the boundary layer greatly affects the loss 
generation. Investigations by Halstead et al. [20] on a 
large-scale LP-turbine use surface mounted hot films to 
acquire detailed information about the quasi-shear 
stress directly on the blade surface. As investigations by 
Cardamone et al. [14] and Schroder [18] indicate, the 
benefit of the wake-boundary layer interaction can be 
used for design procedure of modern gas turbine 
engines with reduced LPT- blade number without 
altering the stage efficiency. 

Most of the above mentioned studies on LP-turbine 
cascade aerodynamics have largely concentrated on the 
measurement of the signals stemming from hot film 
mounted on the suction and pressure surfaces of the 
blades under investigation. Although this technique is 
effective in qualitatively reflecting the interaction of the 
unsteady wake with the boundary layer, because of lack 
of an appropriate calibration method, it is not capable 
of quantifying the surface properties such as the wall 
shear stress. The few boundary layer measurements are 
not comprehensive enough to provide any conclusive 
evidence for interpretation of the boundary layer 
transition and separation processes and their direct 
impact on the profile loss, which is a critical parameter 
for blade design. Furthermore, numerical simulation of 
unsteady LPT-blade aerodynamics using conventional 
turbulence and transition models fails if applied to low 
Reynolds number cases. Recent work presented by 
Cardamone et al. [14] shows that in the steady state 
case at Re= 60,000, the separation is captured, 
however, in the unsteady case, the separation zone is 
not reproduced. 

The objective of the present study, which is the 
first of a series of investigations dealing with specific 
issues of LPT-boundary layer aerodynamics, is to 
provide detailed unsteady boundary flow information to 
understand the underlying physics of the inception, 
onset, and extent of the separation zone. Furthermore, 
the unsteady boundary layer data from the present and 
planned experimental investigations will serve to extend 
the intermittency based unsteady boundary layer 
transition model developed by Schobeiri and his co- 


workers [7, 11, 12] to the LPT-cases, where separation 
occurs on the suction surface at low Reynolds number 
at design and off-design points. Furthermore, the 
experimental results are intended to serve as a 
benchmark data for comparison with numerical 
computation using DNS, LES or RANS-codes. 

It is well known that boundary layer measurement 
is one of the time consuming aerodynamic 
measurements. Any attempt to increase the number of 
parameters to be studied would inevitably result in 
substantial increase of the measurement time. 
Considering this fact, the research facility described in 
[9] and [10] with state-of-the-art instrumentation has 
been substantially modified to study systematically and 
efficiently the influence of periodic unsteady and highly 
turbulent flow on LPT-cascade aerodynamics at the 
design and off-design incidence angles, where Reynolds 
number, wake impingement frequency, free-stream 
turbulence, and the blade solidity can be varied 
independently. 

EXPERIMENTAL RESEARCH FACILITY 

To investigate the effect of unsteady wake flow on 
turbine and compressor cascade aerodynamics, 
particularly on unsteady boundary layer transition, a 
multi-purpose large-scale cascade research facility was 
designed and has been taken into operation since 1993. 
Since the facility in its original configuration is 
described in [9] and [10], only a brief description of the 
modifications and the main components is given below. 
The research facility consists of a large centrifugal air 
supplier, a diffuser, a settling chamber, a nozzle, an 
unsteady wake generator, and a turbine cascade test 
section (Fig. 1). An air supplier with a volumetric flow 
rate of 15 m 3 /s is capable of generating a maximum 
mean velocity of 100 m/s at the test section inlet. The 
settling chamber consists of five screens and one 
honeycomb flow straightener to control the uniformity 
of the flow. 

Two-dimensional periodic unsteady inlet flow is 
simulated by the translational motion of a wake 
generator (see Fig. 1), with a series of cylindrical rods 
attached to two parallel operating timing belts driven by 
an electric motor. To simulate the wake width and 
spacing originating from the trailing edge of rotor 
blades, the diameter and number of rods can be varied. 
The rod diameter, its distance from the LPT-blade 
leading edge, the wake width and the corresponding 
drag coefficient are chosen according to the criteria 
outlined by Schobeiri et al. [21]. 
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(T) Static pressure blade @ Timing belts, rod attachments (g) Inlet nozzle 

(3J Blade with hot film sensors 0. Transition duct (Id) Hydraulic cylinders 

@ Take generating rods (7) Straight duct (11) Pivot point 

@ Wake generator @ Traversing system © Wake generator e-motor 

Fig. 1: Turbine cascade research facility with the components and the adjustable test section. 



Fig. 2: Cascade geometry and stagger angle are listed in Table 1. Number of blades = 5, 
SPB-1 and SPB-2 are blades with static pressure taps, HFB is instrumented with surface 
mounted hot films to be used for future investigations. 


i : i Large silence chamber with 
honeycom and five screens 
(id) Telescope supprt 

(15) Honeycomb flow straightener 
Traversing slots 
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Table 1: Parameters of turbine cascade test section . 


Parameters 

Values 

Parameters 

Values 

Inlet velocity 

V m = 4 m/s 

Inlet turbulence intensity 

Tu in = 1.0% 

Rod translational speed 

U = 5.0 m/s 

Blade Re-number 

Re = 110,000 

Nozzle width 

W = 200.0 mm 

Blade height 

h B = 200 mm 

Blade chord 

c = 203.44 mm 

Blade stagger angle 

o 

'sO 

<N 

II 

>- 

Blade axial chord 

c m = 182.85mm 

Cascade solidity 

o= 1.248 

Blade suction surface length 

L ss = 270.32 mm 

Cascade angle 

6 = 55° 

Cascade flow coefficient 

$ = 0.80 

Cascade spacing 


Inlet air angle. Fig. 2 

a, = 35.0° 

Exit air angle. Fig. 2 

a 2 = 145.0° 

Rod diameter 

D r = 2.0 mm 

Rod distance to lead, edge 

L R = 122 mm 

Cluster 1 (no rod, steady) 

S R = ° ° 

Q - parameter steady case 

o 

© 

II 

Cluster 2 rod spacing 

S K = 160.0 mm 

Q - parameter for cluster 1 

Q= 1.59 

Cluster 3 rod spacing 

S R = 80.0 mm 

Q - parameter for cluster 2 

Q= 3.18 


The belt-pulley system is driven by an electric motor 
and a frequency controller. The wake-passing frequency 
is monitored by a fiber-optic sensor. The sensor also 
serves as the triggering mechanism for data transfer and 
its initialization, which is required for ensemble- 
averaging. This type of wake generator produces clean 
two-dimensional wakes, whose turbulence structure, 
decay and development is to a great extent predictable 
[21]. The facility was used for many unsteady boundary 
layer transition and heat transfer investigations [9-12] 
that serve as benchmark data for validation of 
turbulence and transition models and for general code 
assessments. 

To account for the inlet flow angle and a high flow 
deflection of LPT-cascade, the entire wake generator 
and test section unit including the traversing system 
were modified to allow a precise angle adjustment of 
the cascade relative to the incoming flow. This is done 
by a hydraulic platform, which simultaneously lifts and 
rotates the wake generator and test section unit. The 
unit is then attached to the tunnel exit nozzle with an 
angular accuracy less than 0.05°, which is measured 
electronically. 


The special design of the facility and the length of 
the belts (L^, = 5,000 mm) enables a considerable 
reduction of measurement time. For the present 
investigation, two clusters of rods with constant diameter 
of 2 mm are attached to the belts. The two clusters with 
spacings S R = 160 mm and S R = 80 are separated by a 
section which does not have any rods, simulating steady 
state case (S R = °°). Thus, it is possible to measure 
sequentially the effect of three different spacings at a 
single boundary layer point. To clearly define the 
influence domain of each individual cluster with the 
other one, the clusters are separated by a certain 
distance. Using the triggering system mentioned above 
and a continuous data acquisition, the buffer zones 
between the data clusters are clearly visible. The data 
analysis program removes the buffer zones and extracts 
the data pertaining to each cluster. Comprehensive 
preliminary measurements were c arried out to ensure that 
the data were exactly identical with the data obtained 
with an arrangement of rods filling the entire belt at 
constant spacing. 

The cascade test section, located downstream of the 
wake generator, includes 5 blades with a height of 200.0 
mm and a chord of 203.44 mm, Fig.2. The airfoil used 
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is the Pratt & Whitney “Pak B” airfoil, whose cascade 
geometry is given in Table 1. The cascade test section 
exhibits the essential flow features such as laminar 
boundary layer separation that is inherent to typical 
LPT- blades. The blade geometry was made available 
to NAS Aresearchers and academia to study the specific 
problems of LPT-flow separation, its passive and active 
control and its prevention. As shown in [9], a minimum 
blade number of 5 is necessary and sufficient to secure 
a spatial periodicity for the cascade flow. The 
periodicity was verified by comparing the pressure 
distributions of blade number 2 (SPB-1) and 4 (SPB-2) 
shown in Fig. 2. These blades were specially 
manufactured for measurement of pressure and showed 
identical pressure distributions. A computer controlled 
traversing system is used to measure the inlet velocities 
and turbulence intensities, as well as the boundary 
layers on the suction and pressure surfaces. The 
traversing system is vertically mounted on the 
plexiglass side wall. It consists of a slider and a lead 
screw that is connected to a DC- stepper motor with an 
encoder and decoder. The optical encoder provides a 
continuous feedback to the stepper motor for accurate 
positioning of the probes. The system is capable of 
traversing in small steps up to 2.5 pm, which is 
specifically required for boundary layer investigations 
where the measurement of the laminar sublayer is of 
particular interest. 

INSTRUMENTATION, DATA ACQUISITION, AND 
DATA REDUCTION 

The data acquisition system is controlled by a 
personal computer that includes a 16 channel, 12-bit 
analog-digital f A/D) board (NI, PCI-MIO-16E-1). Time 
dependent velocity signals are obtained by using a 
commercial 3-channel (TSI, IFA-100), constant 
temperature hot-wire anemometer system that has a 
signal conditioner with a variable low pass filter and 
adjustable gain. A Prandtl probe, placed upstream of 
the diffuser, monitors the reference velocity at a fixed 
location. The pneumatic probes are connected to high 
precision differential pressure transducers (MKS 
220CD, range: 0 to lOmmHg) for digital readout. 
Several calibrated thermocouple are placed downstream 
of the test section to constantly monitor the flow 
temperature. The wake generator speed and the passing 
frequency signals of the rods are transmitted by a fiber- 
optic trigger sensor. The passage signals of the rods are 
detected by the sensor using a silver-coated reflective 
paint on one of the belts. This sensor gives an accurate 
readout of the speed of the wake generator and the 
passing frequency of the rods. The signals of the 


pressure transducers, thermocouples, and trigger sensors 
are transmitted to the A/D board and are sampled by the 
computer. The second and fourth blade are each 
instrumented with 48 static pressure taps. The taps are 
connected to a scanivalve, which sequentially transferred 
the pressure signals to one of the transducers that was 
connected to the A/D board. Using Labview software, 
the data acquisition system including the computer and 
the A/D board was upgraded. Two adjacent blades are 
used for boundary layer measurement. 

Steady and unsteady data are taken by angle 
calibrated, custom designed miniature single hot wire 
probes. At each boundary layer position samples were 
taken at a rate of 20 kHz for each of 100 revolutions of 
the wake generator and low pass filtered at 10 kHz. The 
data were ensemble-averaged with respect to the 
rotational period of the wake generator. Before final data 
were taken, the number of samples per revolution and the 
total number of revolutions were varied to determine the 
optimum settings for convergence of the ensemble- 
average. 

For the steady state case, the instantaneous velocity 
components are calculated from the temperature 
compensated instantaneous voltages by using the 
calibration coefficients. The instantaneous velocity can 
be represented in the following form. 

V = V + v (1) 


Where V is the mean (time-averaged) velocity and v is 
the turbulent fluctuation component. The mean velocity, 
also known as the time-average, is given by: 


V 


M 



( 2 ) 


where M is the total number of samples at one boundary 
layer location. The root mean square value of the 
turbulent velocity fluctuation is obtained from the 
instantaneous and mean velocities by: 



and the local turbulence intensity is defined as: 

. , m 

Tu lo = 4*ioo = 4, i~X( v r n 2 *ioo (4) 

V F N M J = 1 

The ensemble-averaged velocity, fluctuation velocity, 
and the turbulence intensity were calculated from the 
instantaneous velocity samples by: 
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m - <m> = ( 5 ) 

Jy j = i 


- <v,(/,)> = 


^E[W - <F A) > ] 2(6) 

iV 7=1 


<v(t)> 

<Tu(t)> ^ <Tu (t)> = xlOO ( 7 ) 

'' ' ' <F(0=V 

where N= 100 is the total number of wake generator 
periods and M the number of samples taken per period. 
<V (t)> ref is the reference ensemble averaged velocity 
for the particular boundary layer traverse. The 
ensemble-averaged boundary layer parameters such as 
displacement thickness <8[>, momentum thickness 
<5 2 >, and shape factor <H 12 > are calculated as follows: 

< 5 > 

5i , <6,> = / (1 - dy (8 ) 

0 e 


<6 > 


8 2 = < 8 2 > = 


/ 


<V> 


(1 


<v> 


H, 


<v> <v> 

e e 

< 8 ,> 


12 


<H„> = 

12 < 8 -> 


) dy (9) 


( 10 ) 


In the above equations the signs < > refer to ensemble 
averaged quantities. For brevity we will drop < > in 
the following discussions. 


EXPERIMENTAL RESULTS AND DISCUSSION 


Detailed surface pressure and boundary layer 
measurements were performed at three different Re- 
numbers. For each Reynolds number three different 
reduced frequencies were applied. For generation of 
unsteady wakes, cylindrical rods with the diameter d R = 
2 mm were chosen to fulfill the similarity criterion that 
requires the generation of a drag coefficient C D that is 
approximately equal to the C D of the turbine blade with 
the chord and spacing given in Table 1 (for details we 
refer to the studies in [21] and [22], 

To accurately account for the unsteadiness caused 
by the frequency of the individual wakes and their 
spacings, the flow velocity, and the cascade parameters, 
we define a reduced frequency Q that includes the 
cascade solidity o, the flow coefficient cp, the blade 
spacing S B , and the rod spacing S R . Many researchers 


have used Strouhal number as the unsteady flow 
parameter, which only includes the speed of the wake 
generator and the inlet velocity. However, the currently 
defined reduced frequency Q is an extension of Strouhal 
number in the sense that it incorporates the rod spacing 
S R and the blade spacing S B in addition to the inlet 
velocity and wake generator speed. For surface pressure 
measurement rods with uniform spacings as specified in 
Table 1 were attached over the entire belt length. For 
boundary layer measurement, however, clusters of rods 
were used, as described previously. 

Surface Pressure Distributions 

Detailed static surface pressure measurements were 
taken at Re=50,000, 75,000, 100,000, and 125,000 using 
the two static pressure blades SPB-1 and SPB-2 shown 
in Fig. 2. The two blades measured identical static 
pressure distributions confirming the spatial periodicity. 
For each Reynolds number three different reduced 
frequencies, namely Q = 0.0, 1.59, and 3.18, are applied 
that correspond to the rod spacings S R = °°, 160 mm, and 
80 mm. The pressure distributions in Fig. 3, which 
pertain to SPB-1, show the results of one steady and two 
unsteady cases. The pressure signals inherently signify 
the time-averaged pressure because of the internal 
pneumatic damping effect of the connecting pipes to the 
transducer. The noticeable deviation in pressure 
distribution between the steady and unsteady cases, 
especially on the suction surface, is due to the drag 
forces caused by the moving rods. The drag forces are 
imposed on the main stream and cause momentum 
deficiency that lead to a reduction of the total and static 
pressure. 

The time-averaged pressure coefficients along the 
pressure and suction surfaces are plotted in Fig. 3. 
Starting with Re=50,000 and steady state with Q=0 
shown in Fig 3a, the suction surface (upper portion), 
exhibits a strong negative pressure gradient. The flow 
accelerates at a relatively steep rate and reaches its 
maximum surface velocity that corresponds to the 
minimum c p =-4.0 at s/s Q = 0.42 . Passing through the 

minimum pressure, the fluid particles within the 
boundary layer encounter a positive pressure gradient 
that causes a sharp deceleration until s/s 0 =0.55 has been 

reached. This point signifies the beginning of the laminar 
boundary layer separation and the onset of a separation 
zone. As seen in the subsequent boundary layer 
discussion, the part of the separation zone characterized 
by a constant c p - plateau extends up to s/s Q =0.746, thus 
occupying more than 19% of the suction surface. 
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Re=50,000 


Re=75,000 




Fig. 3: Static pressure distribution at four different Re-numbers and reduced frequencies Q = 0, 
1,59, 3.18 (no rod, 160 mm, 80 mm), Start= Separation start for steady and unsteady cases, SE= 
Separation end for steady case, UE= Separation end for unsteady cases. 


Passing the plateau, the flow first experiences a second 
sharp deceleration indicative of a process of re- 
attachment followed by a further deceleration at a 
moderate rate. On the pressure surface, the flow 
accelerates at a very slow rate, reaches a minimum 
pressure coefficient auA 0 « 0.42 and continues to 
accelerate until the trailing edge has been reached. 
Unlike the suction surface, the pressure surface 
boundary layer does not encounter any adverse positive 
pressure gradient that triggers separation. However, 
close to the leading edge, a small depression in the 
curve extending from s/s 0 ~ 0.0.8 to 0.24 indicates the 
existence of a small size separation zone that might be 
attributed to a minor inlet flow incidence angle. 

Considering the unsteady case with a reduced 
frequency of Q = 1 .59 that corresponds to a rod spacing 
of S R = 160.0mm, Fig. 3(a) exhibits two distinctive 
characteristics: (1) A noticeable deviation in pressure 
distribution between the steady and unsteady cases. As 
mentioned above, this deviation is attributed to the 
momentum deficiency that leads to a reduction of the 
total and static pressure. (2) For Pe< 125,000, the 


wakes have a substantial reducing impact on the 
streamwise extent of the separation plateau. As seen in 
Fig. 3(a), the trailing edge of the plateau has shifted 
from s/s Q =0.746 to s/s 0 =0.664. This shift has reduced 
the streamwise extent of the separation plateau from 19% 
to less than 1 1% of the suction surface length which is, 
in this particular case, more than 42% reduction in 
streamwise extent of the separation plateau. Although the 
extent of the separation plateau is not necessarily 
identical with the extent of the separation zone, the 
relative reductions presented above adequately reflect 
the relative reduction of the size of the separation zone. 
Doubling the reduced frequency to Q=3.18 by utilizing 
a rod spacing of S R = 80.0 mm causes a slight shift of the 
cp-distribution compared with Q = 1 .59 -case. One should 
bear in mind that pneumatically measured surface 
pressure distribution represents a time integral of the 
pressure events only. Detailed information regarding the 
structure of the separation zone requires a detailed 
unsteady boundary layer or surface pressure 
measurement by fast response probes, as will be 
discussed in the subsequent sections. Increasing the 
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Reynolds number to Re = 75,000. has not brought 
major changes in steady state c -distribution. This is 
also true for the subsequent higher Reynolds number 
cases at steady state. Fig. 3(b,c,d). However, the 
combination of higher Re-number with unsteady wakes 
reveals that the noticeable deviation in pressure 
distribution between the steady and unsteady cases 
discussed above is diminishing with increasing the Re- 
number as shown in Figs. 3(b,c,d). Two counteracting 
factors are contributing to this deviation. The first 
factor is attributed to the momentum deficiency and the 
associated total pressure losses caused by moving 
wakes, as discussed above. The second factor pertains 
to the energizing effect of the impinging wakes on the 
boundary layer. Although the impinging wakes cause 
velocity and momentum deficits, their high turbulence 
intensity vortical cores provide an intensive exchange 
and transfer of mass, momentum, and energy to the 
blade surface, thus energizing the low energetic 
boundary layer. In conj unction with the surface pressure 
distribution, the kinetic energy of the normal velocity 
fluctuation component plays a crucial role. In case of a 
low Re-number flow, the strong damping effect of the 
wall shear stress has the tendency to reduce the normal 
contribution of turbulence kinetic energy, thereby 
diminishing its surface pressure augmenting effect. 
Increasing the Reynolds number results in a decrease of 
the damping effect of the wall shear stress, allowing the 
kinetic energy of the normal velocity fluctuation 
component to increase the surface pressure, thus 
offsetting the wake deficit effects on the pressure 
distribution. This fact is clearly shown in Fig. 
3(a,b,c,d), where the pressure distributions of unsteady 
flow cases at Q = 1.59 and Q=3.18 systematically 
approach the steady state cases ati?e=75,000, 100,000 
and very visibly at Re= 125,000 . It is worth noting that 
the impact of the unsteady wakes on the extent of the 
separation zone is preserved regardless of the Reynolds 
number variation performed in this study. 

Time Averaged Velocity Distributions 

Following the surface pressure investigations that 
mainly addressed the onset and extent of the separation 
zone discussed previously, comprehensive boundary 
layer measurements were performed to identify the 
streamwise and normal extent as well as the 
deformation of the separation zone under unsteady 
wake flow. The steady state case serves as the reference 
configuration. 

Boundary layer profiles were taken for one steady 
and two unsteady inlet flow conditions on the suction 


surface along 31 streamwise locations parallel to the 
cascade front. After completing the velocity 
measurements, the boundary layer coordinates were 
transformed into a blade orthogonal coordinate system. 
Velocities at blade normal positions were obtained by 
interpolating their transformed values. The results 
showed almost no difference between the interpolated 
and non-interpolated velocity data. Experimental 
investigations were performed for three different values 
of Q = 0.0, 1.59, and 3.18. These values cover the 
reduced frequency range encountered in LPT-design and 
off-design operation conditions. 

The effect of wake frequency on time averaged 
velocity and turbulence fluctuation distributions is shown 
in Fig. 4 at four representative streamwise locations for 



Fig. 4: Distribution of time averaged velocity (a) and 
fluctuation rms velocity (b) along the suction surface for 
steady case Q=0 (S R = °°) and unsteady cases Q=1.59 
(S R =160mm) and Q=3.18 (S R =80mm) atRe= 110,000. 
Note the changes in y-scale. 

Re= 110,000. Upstream of the separation zone at 
s/s 0 = 0.52 and its proximity s/sq = 0.588 , the velocity 
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distributions inside the boundary layer are not affected 
by wakes. Inside the separation zone at s/s 0 =0.105 , a 
substantial influence of the wake frequency is observed. 
Higher wake frequency introduces fluctuation kinetic 
energy into the boundary layer which tends to reverse 
the separation tendency. Velocity distributions at 
s/s 0 = 0.767, 0.805, 0.849, and 0.898 (not presented) 
clearly show that the wake impingement shortens the 
streamwise extent of the separation zone compared to 
the steady case. Downstream of the separation zone, 
where the flow is fully reattached, s/s Q = 0.951, the 
impact of wake on the boundary layer is reduced. This 
effect is clearly shown in the velocity distribution at 
s/s Q = 0.951. In accord with the previous 
investigations by Schobeiri et al. [10] on a HP-turbine 
cascade, an increased wake frequency causes 
turbulence fluctuations to rise inside and outside the 
boundary layer as shown in Fig 4(b). However, once the 
boundary layer is re-attached and the velocity 
distribution assumes a turbulent profile, no major 
changes are observed in velocity as well as fluctuation 
distribution. 

Ensemble-Averaged Boundary Layer Velocity 
Distributions 

Figure 5 displays two representative temporal 
ensemble-averaged velocity distributions for (a) steady 
and (b) unsteady flow condition with their characteristic 
features. Both figures show the boundary layer 
development from the freestream to the blade surface at 
a streamwise position of s/s G = 0.0208. Approaching the 
wall surface, both velocities experience a continuous 
deceleration. The velocity gradient in both cases causes 
generation and formation of vortices that transform the 
steady nature of case (a) into an unsteady one as clearly 
demonstrated in Fig. 5(a). The unsteady case displayed 
in Fig 5(b) is characterized by its deterministic temporal 
periodicity. Approaching the wall surface from 
y=10.1mm to 3.45mm, the traveling periodic wake 
experiences a phase shift, while maintaining its 
deterministic nature. However, by penetrating into the 
boundary layer, the interaction between wake and 
boundary layer causes the deterministic nature to 
degenerate into a stochastic one. The results presented 
in Fig. 5 are in full agreement with those discussed in 
[7] and [12], 

Temporal Behavior of the Separation Zone 
Under Unsteady Wake Flow 

Velocity distributions on the suction surface with 
time as the parameter are plotted in Fig. 6. The 



500 550 600 


Time (ms) 

Fig. 5: Ensemble averaged velocity as a function of 
time for (a) steady flow case Q=0 (S R = °°) and (b) 
unsteady case Q=1.59 (S R = 1 60 mm) at s/so=0.0208 
and Re = 110,000. 

nondimensional time (t/r) values are chosen so that they 
represent the temporal states within one full period of 
wake passing. As Figs. 6(a) to 6(e) show, the velocity 
distributions inside and outside the boundary layer at 
fixed s/s 0 -locations experience moderate to pronounced 
changes. Figure 6(a) represents the instantaneous 
velocity distribution upstream of the separation zone 
followed by Figs. 6(b,c,d,e) which represent the velocity 
distributions inside the separation zone. The last figure 
6(f) exhibits the instantaneous velocity distribution 
downstream of the separation zone. In discussing the 
following results, we simultaneously refer to the wake 
distribution as well as the turbulence fluctuation results. 

Figure 6(a) exhibits the velocity distribution on the 
suction surface at s/s„=0.402. At this streamwise 
position, the laminar boundary layer is subjected to a 
strong negative pressure gradient. The boundary layer 
distributions at different (t/r) experience changes in 
magnitude that reflect the corresponding changes of the 
impinging periodic wake velocity. It is worth noting, that 
despite the injection of turbulence kinetic energy by the 
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impinging wakes, no local instantaneous boundary layer 
transition occurs. This is because of the strong negative 
pressure gradient that prevents the boundary layer from 
becoming instantaneously transitional. Instantaneous 
velocity distributions inside the separation zone are 
shown in Figs. 6(b,c,d,e,f). 



V(m/s) V(m/s) V(m/s) 


Fig. 6: Distribution of the ensemble averaged velocity 
development along the suction surface for different s/s„ 
with time t/r as parameter for Q= 1 .59 (S R = 160mm) and 
Re = 1 10,000, note scale change inff). 

As a representative case, we discuss the results plotted 
in Fig. 6(e) at s/s 0 = 0.674. During the time interval 
from tAu close to 0.5 (1.5, 2.5, etc.) to about t/i = 0.75 
(1.75, 2.75 etc.), the separation zone is exposed to the 
wake external flow being under the influence of 
relatively lower turbulence. This flow does not have the 
capability to suppress the separation zone. Thus the 
separation region is clearly shown by the velocity 
distributions at t/r = 0.5 and t/t = 0.75. As the wake 
passes over the blade at s/s 0 = 0.674 introducing high 
turbulence kinetic energy into the boundary layer, the 
boundary layer is energized causing the separation zone 
to partially reduce or disappear. This leads to an 
instantaneous re-attachment. This time interval 
corresponds to the case where the flow is completely 
under the influence of wake and correspondingly the re- 
attached velocity distribution assumes a turbulent 
profile characterized by the curves at t/T =1 .0. t lx = 


0.05, and t h = 0.25 shown in Fig 6(e). To emphasize 
this statement, the steady state velocity distribution at the 
same streamwise position is also plotted in Fig. 6(e) 
using full circles. It shows clearly the separated nature of 
the boundary layer which coincides with the 
instantaneous velocity profile at t/i = 0.5. Intermediate 
times reflect the gradual change between the separation 
and re-attachment as the flow is undergoing the influence 
of the oncoming wake. Moving to the trailing edge of the 
separation zone, at s/s 0 = 0.705, Fig. 6(f), a partial 
reduction in boundary layer thickness as the result of 
wake impingement is visible, however, the separation 
zone does not seem to disappear. 

Temporal-Spatial Resolution of the Separation 
Zone 

To better understand the underlying physics of the 
LPT-flow separation, detailed unsteady flow 
measurements are performed to identify the onset and 
extent of the separation zone discussed previously. The 
separation zone can be thought of as a curve that 
connects the velocity inflection points along the suction 
surface. Starting with a reduced frequency of Q = 1.59 
(S R = 160 mm) atsA 0 =0.520 , Fig. 7(a) exhibits the start 
of the separation zone with a normal extent of less than 
1.0 mm. The impingement of the periodic wake vortical 
core with the high turbulence intensity causes a local 
periodic contraction of the zone in normal direction. 
Convecting downstream, the normal extent increases, 
thus the contraction appears more pronounced. Fig 7 (b, 
c, d). Substantial contraction occurs toward the trailing 
edge of the separation zone as shown in Fig. 7(e, f). 
Similar flow picture is observed when operating at a 
reduced frequency of Q = 3. 18 (S R = 80 mm). 

Studying the temporal distribution of the turbulence 
fluctuations along the streamwise extension of the 
separation zone, it is proposed that, in conjunction with 
the pressure gradient and periodic wakes, there is another 
crucial mechanism responsible for the normal and 
streamwise decrease of the separation zone. This 
mechanism constitutes a combination of the high 
turbulence fluctuation level and its gradient. It is the 
temporal gradient of the turbulence fluctuation, or more 
precisely, the fluctuation acceleration dv rms /dt that 
provides higher momentum and energy transfer into the 
boundary layer energizing the separation zone and 
causing it to partially or entirely disappear. 
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y(mm) y(mm) Y(mm) 





S R =160mm, s/s o =0.617 (d) 


V(m/s) 

0 8.14 
7.79 
7.43 
7.08 
6.72 
6.36 
6.01 
5.65 
5.30 
4.94 
4.58 
— | 423 
3.87 
3.52 
3.16 
2.00 


V(nVs) 

0 8.01 
7.64 
7.28 
6.92 
6.55 
3 6.19 
5.83 
5.46 
5.10 
4.74 
4.38 
4.01 
3.65 
3.29 
2.92 
2.00 


S R =160mm, s/s o =0.588 



Fig. 7: Contour plot of the ensemble averaged velocity distribution showing the effect of periodic wakes on 
the separation zone at different streamwise positions and Re = 1 10,000. 
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For better understanding this phenomenon, we present 
the wake and the fluctuation velocity in Fig. 8 for the 
streamwise position .v/.s- () = 0.651 indicated in Fig. 7(e). 
For the sake of clarity, we scale up the fluctuation 
velocity with a factor of 4 and choose an arbitrary 
normal position of y = 2.85mm to be sufficiently above 
the separation zone. 



Fig. 8: Wake velocity and fluctuation rms distribution, 
Re= 110,000. 


fluctuation within the vortical core starts to increase while 
the velocity continuously decreases. This process 
continues until the end contraction at //x= 1.41 has been 
reached. Thereafter, the separation zone is subjected to a 
process of intensive exchange of momentum and energy 
that causes the separation to diminish, as shown in Fig. 9. 
The process of separation contraction, suppression, and 
regeneration is summarized in Fig. 10. It shows more 



Fig. 10: Separation contraction, suppression, and 
regeneration, Re=l 10,000. 


Figure 8 exhibits two distinct regions: (1) a wake 
vortical core, occupied by vortices that originate from 
the moving cylindrical rods and generate high 
turbulence fluctuations, and (2) a wake external region 
between the adjacent vortical cores with relatively low 
turbulence activities. The wake configuration is 
asymmetric as discussed in [21. Figure 7(e) is enlarged 
in Fig. 9 to reveal further details. 



Fig. 9: Separation zone, definition of contraction 
begin, contraction end, suppression, and regeneration. 


As Fig. 9 shows, the separation zone starts to contract 
at f/x=1.25 (2.25 etc.) This point coincides with the 
streamwise position of the velocity maximum, which 
exactly corresponds to the position of the fluctuation 
minimum, as shown in Figs. 8 and 10. At this point, the 


details of separation contraction and suppression. In this 
context, it is necessary to subdivide the vortical core 
shown in Fig. 8 into four distinct regions, separated by 
thick dashed lines, as presented in Fig. 10. Region (a) is 
characterized by the initial positive gradient of the 
fluctuation dv/dt>0 marked with an upward arrow. 
Region (b) represents the substantial part of the vortical 
core with an intense turbulence activity. Region (c) serves 
as a transition region between region (b) and the relatively 
calm region (d) characterized by dv/dt< 0 . 

For an initial fluctuation gradient dv/dt> 0, the 
separation zone begins to contract at f/x= 1.25 . This initial 
gradient is crucial for initiating the contraction process. 
Once the maximum fluctuation velocity with the temporal 
gradient dv/dt = 0 at f/x= 1.41 is reached (region a), the 
process of energizing the separation zone starts 
transferring momentum and energy into the separation 
zone, thereby preventing its regeneration (region b). 
Passing the transition region (c), the process of 
suppression continues until the end of the vortical region 
at f/x=2.0 is reached. At this point the external wake 
region with its low turbulence content arrives causing a 
regeneration of the separation zone, thus reversing the 
entire suppression process. 

While turbulence fluctuation expressed in terms of 
higher turbulence intensity is well known for influencing 
the flow separation, its gradient enhances the effect of 
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delaying the onset and reducing the extent of the 
separation zone. The fluctuation gradient is an inherent 
feature of the incoming periodic wake flow and does 
not exist in a statistically steady flow that might have 
high turbulence intensity. 

The results clearly indicate that for the particular 
blade under investigation, one has to deal with a large 
separation on the suction surface. These observations in 
comparison with the steady state reference case suggest 
that, once a separation zone is manifested, its size can 
be significantly reduced by periodic wake impingement, 
but it cannot be completely removed. The results 
presented here are valid for blades with a similar 
pressure distribution discussed earlier. Since the onset 
and extent of the separation zone is uniquely associated 
with the pressure gradient, blades can be designed with 



Fig. 11: Ensemble-averaged relative momentum 
thickness distribution along the suction surface for 
Q =1.59 ( S R = 1 60 mm) and Re = 110,000. 


less local adverse pressure gradient, whose separation 
onset can completely be suppressed by impinging wakes. 

Boundary Layer Ensemble Averaged Integral 
Quantities 

The integral parameters, such as momentum thickness 
and shape factor, are of particular interest to a turbine 
designer, since they provide an accurate first estimation of 
the quality of the designed blade. The ensemble-averaged 
distributions of the momentum deficiency thickness and 
shape factor for the suction surface are shown in Fig. 
ll(a,b,c,) for Q = 1.59 at different s/s 0 -locations . The 
momentum thickness values are non-dimensionalized with 
respect to the value corresponding to the steady case with 
Q = 0. The period t represents the wake -passing period 
that is specific to the individual wake generating cluster, 
which is characterized by the Q - value under 



Fig. 12: Ensemble averaged relative form parameter 
distribution along the suction surface for Q = 1.59 
(S R = 160mm), Re = 110,000. 
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investigation. The periodic behavior of the ensemble- 
averaged momentum thickness over the entire suction 
surface as a result of the embedded periodic wake flow 
is clearly visible from Fig. 11. It represents the 
momentum thickness behavior at different s/s 0 - 
locations upstream, within and outside the separation 
zone. The relative momentum thickness distribution 
upstream of the separation zone in Fig. 1 1(a) integrally 
exhibits a slight increase, whereas inside the zone 
shown in Fig. 1 1(b), a substantial decrease is apparent. 
At the immediate vicinity of the separation zone trailing 
edge, s/s 0 = 0.705 close tore-attachment. Fig. 1 1(c), the 
momentum thickness experiences a noticeable increase, 
which by convecting downstream decreases again and 
approaches the integral values that are close to the 
steady state case. 

The ensemble averaged relative shape factor H 12 - 
distributions on the suction surface at the same 
streamwise locations are plotted in Fig. 12(a,b,c). 
Upstream of the separation zone. Fig 12(a) they 
experience a similar periodic change with an average 
value that is close to the value of the steady case. 
Moving into the separation zone. Fig. 12(b,c), each 
streamwise location presents its own form parameter 
that is specific to the velocity distribution we discussed. 

CONCLUSIONS 

A detailed experimental study on the behavior of 
the separation zone on the suction surface of a highly 
loaded LPT-blade under periodic unsteady wake flow 
was presented. One steady and two different unsteady 
inlet wake flow conditions with the corresponding 
passing frequencies, wake velocity and turbulence 
intensities were investigated utilizing a new large-scale, 
subsonic research facility. The results of the unsteady 
boundary layer measurements were presented in 
ensemble-averaged, and contour plot forms. Surface 
pressure measurements were performed at Re= 50,000, 
75,000, 100,000, and 125,000. At each Reynolds 
number, one steady and two periodic unsteady inlet 
flow measurements were performed. The surface 
pressure distribution showed no major changes with 
respect to the above Re -number changes. Noticeable 
changes occurred while operating at unsteady flow 
conditions. Detailed unsteady boundary layer 
measurement identified the onset and extent of the 
separation zone as well as its behavior under unsteady 
wake flow. Passing the wake flow with its highly 
turbulent vortical core over the separation region 
caused a periodic contraction and expansion of the 
separation zone. It was proposed that, in conjunction 


with the pressure gradient and periodic wakes, the 
temporal gradient of the turbulence fluctuation, or more 
precisely the fluctuation acceleration dv r /dt provides 
higher momentum and energy transfer into the boundary 
layer energizing the separation zone and causing it to 
partially or entirely disappear. We found that 
fordv /dt> 0 , the separation zone starts to contract 
whereas for OV rms /dt < 0 it gradually assumes the shape 
before the contraction. The existence of higher turbulence 
fluctuations expressed in terms of higher turbulence 
intensity is well known for influencing the flow 
separation; its gradient is of crucial importance in 
suppressing or preventing the onset and the extent of the 
separation zone. The fluctuation gradient is an inherent 
feature of the incoming periodic wake flow and does not 
exist in a statistically steady flow that might have high 
turbulence intensity. 

UNCERTAINTY ANALYSIS 

The Kline and McClintock [23] uncertainty analysis 
method was used to determine the uncertainty in the 
velocity after calibration and data reduction for the single- 
wire probe. The Kline and McClintock method determines 
the uncertainty for a 95% confidence level. The 

uncertainty in velocity for the single-wire probe after data 
reduction is given in Table 2. As shown, the uncertainty 
in velocity increases as the flow velocity decreases. This 
is due to the pneumatic pressure transducer having a large 
uncertainty during calibration. 


Table 2 : U ncertainty in velocity measurement for hot-wire 
probe. 


U (m/s) 

3 

5 

12 

wU/U (%)„,, 

5.78 

2.41 

1.40 
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